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Abstract 

Both the robust INTEGRAL 511 keV gamma-ray line and the recent tentative hint of the 135 
GeV gamma-ray hne from Fermi-LAT have similar signal morphologies, and may be produced from 
the same dark matter annihilation. Motivated by this observation, we construct a dark matter 
model to explain both signals and to accommodate the two required annihilation cross sections 
that are different by more than six orders of magnitude. In our model, to generate the low-energy 
positrons for INTEGRAL, dark matter particles annihilate into a complex scalar that couples to 
photon via a charge-radius operator. The complex scalar contains an excited state decaying into 
the ground state plus an off-shell photon to generate a pair of positron and electron. Two charged 
particles with non-degenerate masses are necessary for generating this charge-radius operator. One 
charged particle is predicted to be long-lived and have a mass around 3.8 TeV to explain the dark 
matter thermal relic abundance from its late decay. The other charged particle is predicted to 
have a mass below 1 TeV given the ratio of the two signal cross sections. The 14 TeV LHC will 
concretely test the main parameter space of this lighter charged particle. 



1 Introduction 



Although dark matter serves as the dominant component of matter in our universe, its various prop- 
erties remain unknown. From astrophysical evidence, there is no doubt that dark matter can interact 
with the Standard Model (SM) particles through gravitational interaction. However, whether there 
are additional interactions between dark matter and SM particles is still a mystery to us. Among 
several approaches to search for dark matter particles, measuring the cosmic ray spectrum provides 
the indirect detection of dark matter. Observing a high-energy gamma-ray line has long been believed 
to be the "smoking gun" of the dark matter detection [THS]. Furthermore, the propagation of ener- 
getic photons in our Galaxy is less affected by the interstellar gas or Galactic magnetic field. The 
gamma-ray line signal can even provide the dark matter geometrical profile in our Galaxy. 

The detection of celestial gamma-ray line at 511 keV from the inner galaxy, which is believed 
to be caused by e~^e~ annihilations, was first reported by [7J and later confirmed by [HHllj. The 
total flux of the 511 keV line has been estimated to be around 2 x 10^'^ cm^^ s^^ [12]. About 97% 
e~^e~ annihilations proceed through the intermediate state of a positronium atom, and 25% of these 
annihilations with opposite spin of and e~ can produce 511 keV line emission |13lll4j . Although 
this gamma-ray line has been known for decades, the identification of the positron source remains 
undetermined. Different astrophysical sources have been suggested during the years, but each of the 
models faces various challenge to explain the observations consistently. The relatively high ratio of the 
bulge to disk 511 keV emission towards the inner Galaxy seems against its origin from hypernovae and 
gamma ray bursts, while the constraints on the production rate of high energy positrons also disfavors 
millisecond pulsars, as well as proton-proton collisions from e.g. microquasars, low luminosity X-ray 
binary jets, and the central supermassive black hole. Furthermore, pulsars, magnetars, and Galactic 
cosmic rays are not favored as major sources to the observed 511 keV from the bulge, and stringent 
constraints on these origin of the 511 keV line was suggested |15pi6j. 

Besides these astrophysical suggestions, the possibility that dark matter may create the 511 keV 
gamma-ray line has been widely discussed, mainly motivated by the rather spheroidal, symmetric, 
bulge-centered morphology. The lack of higher energy gamma ray requires the injection energy of 
positrons to be below ~ 3 MeV [15]. This motivates studies for both MeV-scale dark matter mod- 
els |17H22j and TeV-scale dark matter models with a MeV mass splitting among different dark matter 
states |23H26j . Since TeV-scale dark matter with electroweak interaction strength can naturally gives 
correct thermal relic abundance, those models are more favored. Interestingly, the morphology of the 
511 keV signal profile has a peaked structure around the Galactic center, and the sharpness of the 



1 



peak prefers to have dark matter annihilation rather than decaying as an explanation [T2]. Thus we 
focus on the heavy dark matter scenario, and try to explain the 511 keV INTEGRAL signal via dark 
matter annihilation. 

One popular dark matter model to explain the INTEGRAL signal is the excited dark matter 
model with an MeV mass splitting [21]. This class of models suffer from the requirement of a large 
kinematic energy of dark matter to excite the ground state, hence relying on the Boltzmann tail of dark 
matter velocity distribution. It is under a debate whether the excited dark matter model can generate 
enough positrons to explain the large gamma-ray flux for the INTEGRAL data |27y28j. For the 100 
GeV dark matter mass region that we will consider in this paper, the situation is worse, because it 
requires a higher velocity to obtain enough kinematic energy comparing to a TeV mass dark matter. 
In our paper, we will address this problem and propose a new scenario of the Down-scattering excited 
Dark Matter (DeDM) to solve the Boltzmann suppression problem of the vanilla excited dark matter 
models. 

More recently, the hint for another gamma-ray line around 130 GeV from the Galactic center 
has been suggested by analyzing the public data from Fermi Gamma-ray Space Telescope (Fermi- 
LAT) j29tl30]. The hint becomes even stronger with the template fitting approach, which takes into 
account the spatial distribution of the LAT events towards the inner Galaxy along with the spectral 
information [31]. Fermi-LAT Collaboration has confirmed the hint of the peak at ~ 130 GeV using 
Pass 7 data. The peak shifts to a higher mass at ~ 135 GeV and the significance becomes weaker 
using the reprocessed Pass 7 data [32] . Such high energy gamma-ray line emission has been considered 
as a clean signature from dark matter annihilations. Many dark matter models have been constructed 
to explain the 130(135) GeV gamma-ray line feature (see [S^ and references therein). 

The morphology of the INTEGRAL 511 keV and Fermi-LAT 130(135) GeV line shares similar 
structures: (1) the signal events concentrate at the center of the Galaxy with non-disk like dis- 
tributions; (2) after smoothing Fermi-LAT signal with INTEGRAL'S angular resolution, they have 
comparable full widths at half maximum (FWHM) in both the longitudinal and latitude directions. 
This motivates us to explain both signals using same dark matter particle in our universe. The fittings 
for both signals prefer annihilation rather than decaying |12ll34j . Having worked out the required an- 
nihilation cross sections, we find that the INTEGRAL 511 keV cross section is six to seven orders of 
magnitude larger than that of the Fermi-LAT 130 GeV line. This large hierarchy of cross sections sets 
a challenge when constructing a detailed model. However, the order of magnitude is comparable with 
an electromagnetic loop factor of ©(a^/vr^) if the INTEGRAL and Fermi-LAT signals are coming 
from tree-level and loop-level processes, respectively. This serves as a clue for our model building. 
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Our paper is organized as following. In ^ we emphasize the similarities of morphologies for both 
signals and work out the required cross sections. In f|3l we propose our model, Down-scattering excited 
Dark Matter model, to explain both signals. In section ^3.11 we first provide a general operator analysis 
to illustrate the essence of our model and calculate the scales of cutoffs of the effective operators. Then 
we build up a concrete UV-completion for the operator analysis in §3.21 In §3.3^ we discuss the thermal 
history of our model. One charged particle needs to be long-lived in our UV model, so that we have 
a semi-natural model to explain the final dark matter relic abundance. 

2 Experimental Data 

In this section, we discuss the INTEGRAL and Fermi-LAT oberservations in more detail. Since 
photon is not much affected during its propagation in the Galaxy, photon coming from dark matter 
annihilation can be used to determine the dark matter distribution in our Galaxy. However, there are 
subtleties on how to map the INTEGRAL 511 keV gamma-ray line signal profile to the dark matter 
distribution profile. This is because the low-energy positrons that are generated from dark matter 
particles can propagate through the interstellar medium and annihilate with electrons to photons 
away from the production site and bias the inferred dark matter distribution from the 511 keV line 
morphology. In this paper, we assume that the positron propagation is negligible comparing to the 
spatial resolution of the INTEGRAL, thus the dark matter profile can be estimated by measuring 511 
keV emission morphology. On the other hand, the Fermi-LAT 130(135) GeV photons could directly 
be generated from dark matter particles, and its morphology can therefore tell us the dark matter 
profile. 

We first compare the morphologies of the INTEGRAL 511 keV and Fermi-LAT 130(135) GeV 
lines. After smoothing the Fermi-LAT 130(135) GeV line using the angular resolution of INTEGRAL, 
we find the spatial distributions are comparable to each other. Furthermore, assuming that both 
signals are generated by dark matter annihilation, we estimate the annihilation cross sections for the 
two processes. They will serve as inputs for our model building in the rest of the paper. 

2.1 Experimental Data 

Thanks to a coded mask located 1.7 m above the detector plane and a specific dithering observational 
strategy, the spectrometer (SPI) onboard the INTEGRAL observatory can image the sky with a 
spatial resolution of ~ 2.6° (FWHM). Based on observations recorded from February 22nd, 2003 to 
January, 2nd 2009, the study in [35j has obtained the morphology of the 511 keV line towards the 
inner Galaxy. In Figure [U we compare the intensity of the 511 keV gamma-ray line as a function of 
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Figure 1: Comparison of the INTEGRAL 511 keV line profile [35] and Fermi-LAT 130 GeV line 
profile [36] from the Galacitc center, on longitudinal (left) and latitudinal (right) projections. The 
black dashed line shows the 511 keV line profile measured by INTEGRAL, and the black solid line 
shows the fitted Gaussian. The shaded histogram shows the 130 GeV line profile from 3.7 years of 
Fermi-LAT data. The red solid line shows the best fitted Gaussian of the 130 GeV line, which is the 
same as Fig. 15 in [31]. The green solid line illustrates the 130 GeV line profile smoothed by SPI 2.6° 
FWHM beam. For the INTEGRAL data, the vertical axis has an arbitrary scale in this plot. 



Galactic longitude and latitude with the 130 GeV line profile obtained by fitting 3.7 years Fermi-LAT 
observations [36]. Especially, the dark green line shows the 130 GeV line profile further smoothed by 
SPI 2.6° FWHM beam. 

Interestingly, both longitudinal and latitude distributions of INTEGRAL are comparable to those 
of Fermi-LAT after smoothing. Furthermore, both distributions show the tendency to be off-center 
in the negative longitudinal direction [3HI35j. These similarities motivate the attempts to build one 
dark matter model to explain both these two signals. 

2.2 Dark matter annihilation cross sections for INTEGRAL and Fermi-LAT 

As discussed in the previous section, both the 511 keV line and the Fermi-LAT 130(135) GeV line 
could potentially be explained by dark matter annihilation. In this section, we estimate the required 
annihilation cross sections for both experiments. 

The gamma-ray line intensity in a given direction provided by dark matter annihilation is the 
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line-of-sight integral of the squared dark matter number density along that given direction 

jE-^-dx Neon. ^ 6iE - E,) , (1) 



with the J- factor defined as: 



J= I db I d£ I ^ cos6 



p(r) 
PQ 



(2) 



Rq 

where / and b are longitude and latitude, and the integral of s is along the line of sight. Here, 
Rq ~ 8.5 kpc is the distance from the Sun to the galactic center; p{r) is the dark matter halo profile; 
Pq ~ 0.4 GeV cm~'^ is the often-used dark matter density in the Solar system [37]; the relation 
between r and s is = + Rq — 2sRq cos £ cos 6; Ny {N^+) is the number of photons (positrons) 
generated from each dark matter annihilation hard process; is the dark matter mass; {avr)e+ and 
{crvr)^ are the annihilation cross sections. We define = 1 for self-conjugated dark matter, e.g. a 
real scalar or a Majorana fermion, and d^ = ^ for a complex scalar or a Dirac fermion. A'^conv is the 
number of monochromatic photons that the final states could convert to. For Fermi-LAT, Nconv = !> 
since we assume that only monochromatic photons are produced in the hard process. For INTEGRAL, 
observations suggest that about 97% of positrons annihilate through positronium formation [38 1. Only 
1/4 of annihilation takes place in the parapositronium state, which produces two 511 keV photons. 
So, we have A^conv ~ 0.55 for INTEGRAL. 

We consider both the Einasto and the Navarro- Frenk- White (NFW) dark matter profile 

-I 3-a 

(3) 



-l\(j-Y-(!:eY] /rQ\c^ 



1 + rQ/rs 



1 + r/rs 

with = 20 kpc and a = 0.17 for Einasto [39] and q = 1 for NFW [IQ]. Using the fitted fluxes for the 
INTEGRAL signal (the dark matter -|-disk ones) in Ref. [H], we obtain the annihilation cross sections 
as0 

(^Tt;r)^,511,Em(NFW) = X 1.5(0.28) X 10^ X pb • C . (4) 

For the Fermi-LAT 130(135) GeV gamma line, we use the fitted fiuxes from Ref. [29] for both profiles 
to calculate the cross sections, 

(<^^^r)7,135,Em(NFW) = X 0.42(0.76) X W'^ X ( ^qq g^y j Pb " C . (5) 



^Here, we use different parameters for dark matter profiles compared to the ones in Ref. |41] . We simply rescale their 
signal flux by the ratio of J functions, which could bring an uncertainty of 0(1). 
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To quantify the ratio of the required cross sections for two experimental results, we define R, 
((Tt>,.)^^i35/((TWr)'y,5ii. Taking = = 2, we have the experimentally measured ratios as 



where clearly show a large hierarchy for the two required cross sections. We want to also stress that the 
astrophysical uncertainties are fairly large and a global fit by combining the INTEGRAL and Fermi- 
LAT might bring the uncertainty down. The cross section ratio between these two expertiments is 
0(10^^ ~ 10^''). This will be the input for model building in latter sections. Interestingly, this ratio 
is comparable to the square of the electromagnetic loop factor (a^/vr^) ~ 6 x 10~^, which implies these 
two experimental results may be related by a loop with two electromagnetic vertices. It serves as a 
clue for model building. 

3 Down-scattering excited Dark Matter 

There are several interesting features required to construct dark matter models if both signatures are 
to be explained by the same dark matter particle with a mass at the 100 GeV scale. 

• The required cross section for the INTEGRAL data is amazingly large. For a simplest estimation 
on the annihilation rate, we get av ~ l/(47rm^) ~ 3 x 10^ pb • c for around 100 GeV. This 
estimation is three orders of magnitude smaller than the required one. Additional mechanisms 
are therefore required to increase the annihilation rate. There are several ways to achieve this 
and we pay special attention on the resonance enhancement |42H47j . 

• To explain the INTEGRAL data, primary positron injections from dark matter are required. 
Since we don't see any excess for other cosmic rays, the underlying dark matter model should 
be arranged to treat positron/electron differently from other particles. In principle, this can be 
achieved either from kinematic constraints or symmetry reasons. 

• The ratio of the two cross sections is {(TVr)is5/ (crvr) 511 10^'' or 10~^. The dark matter model 
should also provide a natural explanation for this hierarchy of two cross sections. 

• The model should provide correct amount of dark matter relic abundance to be consistent with 
observation. 

In the following, we will provide a particle physics model to incorporate all above four ingredients. 
Specifically, we will use a resonance particle in the s-channel to increase the dark matter annihilation 



(i?5n)exp,Ein « 2.9 X 10 



(-^51l)exp,NFW ~ 2.7 X 10 
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cross section required for 511 keV gamma-ray line. The kinematic constraints from a small mass split- 
ting will be used in this paper to select positron/electron as the signals from dark matter annihilation. 
Instead of introducing a light mediator, e.g. dark photon, for the dark matter sector connecting to 
the positron/electron, we use photon as a more natural mediator to achieve this goal. Noticing that 
a neutral scalar field cannot decay into another neutral scalar field plus one on-shell photon, which is 
the reason why T[(n -|- 1)5] T[nS] -|- 7, a neutral scalar coupling to photon with the charge-radius 
operator can naturally induce without generating a photon signal in the meanwhile. The mass 
difference of the two scalar fields is chosen to be small such that the kinematic energy of is small 
enough to be consistent with observation. To explain the ratio of the cross sections, we will have the 
cross section for INTEGRAL to be controlled by coefficients of renormalizable operators, while loop- 
generated higher-dimensional operators for Fermi-LAT. We will first perform an operator analysis and 
then provide a UV-complete model. 

3.1 Operator analysis 

We introduce one Dirac fermion x oii^ complex scalar field $ = {(f)i +i(/>2)/\/2 in the dark matter 
sector. Both x ^-iid 4'i ^-re stable particles and coexist in our current Universe. In our study, we will 
assume that the dark matter component x occupies the majority of the dark matter energy, but we 
will come back to discuss the relative relic abundances of them later. The interactions of the dark 
matter sector to the SM particles are described by the following set of effective operators 

- £ D a, X7^X 5 + 5 + ^ 5 F^^F^'^ + d.^d^^^F^^ , (7) 

where we implicitly assume that the higher-dimensional operators can be generated at one-loop level. 
The annihilation of x's is through exchanging the real scalar S in the s-channel. For the INTEGRAL 
data, a small mass scale at around 1 MeV is required to generate positrons almost at rest. In our 
model, we introduce this small mass scale as the mass splitting of 4>i and 02 from # = ((/>! -|- (j)2)/V^ 
such that 5 = m^p^ — fn(j)^ <C m^^ and J ~ 1 MeV. Noticing that the parameter 5 explicitly breaks the 
global U{l)ci), so the smallness of 5 is technically natural. Expanding the last operator in terms of ^\ 
and 02, we have 

.^^,id,4>M.F^'^ . (8) 

Using the equation of motion, one can rewrite the above operator as 02^1/01^//-^'*'^ = 02^i^0ij'^- This 
indicates that 02 cannot decay to a mass-on-shell photon. For 2me < 5 < 2m^, we have the leading 
decay channel of 02 as 

02 ^ 01 + 7* ^ 01 + e+ e" . (9) 
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Figure 2: The Feynman diagrams for INTEGRAL (left) and Fermi-LAT (right). 



Photon, naturally, behaves as a mediator for the dark matter sector to generate positrons. 

The processes to generate positrons for INTEGRAL and photons for Fermi-LAT are shown in 
Fig. O where the solid thick points indicate higher-dimensional operators for those vertices. Although 
it looks like that the relative cross sections for those two processes are unrelated to each other, we 
will show in a concrete renormalizable model that the overall cross sections could have a relation in 
^3.2[ In order to generate slowly moving positron from dark matter annihilation, as preferred from 
the INTEGRAL data, there are two conditions required: (1) the mass splitting 5 should be close to 
2me; (2) (j)2 cannot have a large boost. The first condition can be satisfied by choosing 5 > 2me. The 
second condition can be arranged by choosing m^^ < m^. 

We first calculate the annihilation cross section for INTEGRAL. Using the interactions of 02 in 
Eq. d?]), one gets the annihilation cross section of XX ~^ 4'2(t>2 at leading order in v^. as 



{'^^r)^2<f>2 (4^2 + ^2^2 _ ^2 )2 + ^2 ^2 ^ ^ ^2' ' (^0) 

We are interested in the parameter space with 2my_ > ms > 2m02 • The decay width of S is calculated 
to be Ts^2Tp +rl with 



/2 Ami 



= _J1 , 1 ^ (11) 

~ 64 7r3M2- ^^^^ 



Here we treat the decay width of 5 — t- 4)i(j)i to be approximately the same as F 



<i>2 
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Figure 3: Left panel: the contour plots of the annihilation rates in pb-c for and the mass difference 
(2m^ — ms). The red and dotted line is the required cross section to explain the INTEGRAL data for 
the Einasto profile, while the red and dotdashed line is for the NEW profile in Eq. Right panel: 
the same as the left one but in terms of m^^ and {2m^ — ms)- 



Eor INTEGRAL, we need to calculate the velocity-averaged annihilation rate, which is given by 

{aVr)511 = {{(^Vr),l}2(t>2M) = —?= / VrdVre~'''"'-^^{aVr),j)2<l)2M , (13) 

v47r J 

where x = rn^/T = v^'^ with vo determining the variance of the Gaussian dark matter velocity 
distribution. In our numerical calculation, we neglect the upper limit of the integration, which is 
controlled by the escaping velocity of dark matter in the galaxy and has only a small effect on our 
final results. 

In Fig. [31 we show the contours of the annihilation rates of {crvr)5ii in terms of /i and C^rriy. — ms) 
in the left panel, also and {2m^ — ms) in the right panel. To obtain a large annihilation rate 
around 10^ pb to explain the INTEGRAL data, the resonance mass has to be very close to twice of 
the dark matter mass. The mass splitting should be a few GeV for the parameter /u ~ 100 GeV. In 
Fig. El we only presented the results for ms ^ ^m^^,. The case with ms > 2m^ has an additional 
contribution to the total width of S from S — t- xXj similar results. From the right panel of 

Fig.[3l we can see that the annihilation rate is insensitive to m^f^^ except for the region with m^^ ~ m^. 
One might think that 02 can be as light as possible. However, a light (j)2 generated from dark matter 
annihilation can have a large Lorentz boost. As a consequence, e+ from <f>2 decays is also boosted and 
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Figure 4: The ratio of annihilation rates for the Fermi-LAT and INTEGRAL signals as a function of 
the cutoff of the higher-dimensional operators. Here, we choose = 100 GeV. The dotted and dot- 
dashed lines (red) indicate the approximate value from experimental measurements for two different 
dark matter profiles. 

too energetic to explain the INTEGRAL data [12]. Therefore, we restrict the parameter space in our 
later study to have m^^ at least above 50 GeV. 

For the Fermi-LAT signal, instead of obtaining the absolute annihilation rate, we calculate the ratio 
of the required signal strengths for Fermi-LAT and INTEGRAL. The ratio is equal to the branching 
ratio of the two decay channels of S in Eq. ()12p . assuming that the additional contribution from the 
process — >■ 77 is small. By taking the ratio, the dependence on the resonance propagator is 
cancelled and we have 



We show this ratio of the annihilation rates in Fig. S] by fixing = 100 GeV. From Fig. [U we can 
find that the cutoff of the operators in Eq. ([7|) should be several TeV. 

From the effective operator analysis in this section, we have seen that it is possible to explain 
the required annihilation rates for both INTEGRAL and Fermi-LAT. Our model is economical in a 
sense that only a few operators and a small number of degrees of freedom are required to explain the 
data. On the other hand, we should also admit that the resonance requirement of C^m^ — ms) ^ ms 
is a tuning point of the parameter space of the current model. Additional ingredients are therefore 
required to explain this delicate mass relation. We leave this direction of exploration to a future study. 
Here we emphasize that the ratio of INTEGRAL and Fermi-LAT signals are independent on the way 
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Table 1: Matter content and corresponding charge assignments. The global symmetry U{1)^ is only 
an approximate one. The small mass splitting of (pi and 02 breaks it. 

we enhance the annihilation cross sections. Thus one can attach the rest of the model to any other 
ways of enhancement, e.g. a light mediator in the t-channel plus the Sommerfeld enhancement |48j . 
In the next section, we construct a renormalizable model to UV complete the Lagrangian in Eq. ([7]) 
and explain the common origin of the last two operators in Eq. ([7|). 

3.2 Renormalizable model 

One way to UV complete the effective Lagrangian in the previous section is to introduce electro- 
magnetic charged states to connect the dark matter sector to photon. In order to have the state 
(pi stable, at least two charged particles are required to preserve the discrete symmetry associated 
with (pi. As one example, we introduce two charged complex scalar fields, Xi and X2. One could 
also study fermionic charged states in the same procedure. Under U{1)y or U{l)em after electroweak 
symmetry breaking, Xi and X2 have charge one. The global symmetries that we introduce contain 
a Z2 symmetry responsible for the stability of the dark matter particles and a f/(l)</) protecting the 
mass degeneracy of (pi and (p2- We show the field content and symmetries in Tabled) 

Based on the symmetries in Table [U we have the following subset of operators allowed by the 
symmetries, 

-CD iX^xi^xS + fiS^^<i> + ni{<i>XiXl + <i>^xlX2) + H2SX1XI + fisSX2Xl +Xim^HW 
+ 1^2 ^2 ^ ^2^^t + ij„^5($2^^t2) + mj^^XiXl + mj^^X2Xl 

+ X'^Xixeii + h.c.. (15) 

Here we only list the operators which are relevant to the processes we concern in this paper. Especially 
the operator SHH^ is neglected, which is assumed to have a small coefficient. The last operator 
generically introduces lepton flavor violation processes, so the couplings A* (i is the flavor index) 
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should be small. 

We first note that the vertices {^XiX^ + ^^xlX2) could generate the charge radius operator for 
the $ field as shown in the last operator in Eq. ([7|). After a calculation of the triangle diagram with 
Xi and X2 propagating in the loop, we get the Feynman rule of the following operator d^^d,^^'^ F'^'^ = 



^TT^ dx dy 
32 TT^ Jq Jo 



(1 _ 2x) fcf + (1 - 2y) fcg 



— X — y)mj^^ + {x + y) mj^^ + {xki +y ^2)^ — xkf — y k2 



m 



where ki and k2 are momenta of 0i and (f)2 with opposite directions towards the vertex. In the limit 
^> kf,k2, we can match to the coefficient of the effective operator in Eq. ([7]) as 

A$_ ^ ^1 - m%^) - 2{m\^ + 4m^^m^^ + m^J log jmxjmx^)] 

We notice that the above formula vanishes when mxi = rnx2- This can be understood by the 
enhanced discrete symmetry, <I> — )• <1>t,Xi -h- X2, in the Lagrangian when Xi and X2 have degenerate 
massesl^ The charge-radius operator violates this discrete symmetry, thus cannot be generated when 
'^Xi = ■ Another more intuitive explanation is to think <I> as a composite particle of X^ and X2 ■ 
If the mass of X2 is much heavier than X^ , one can treat X^ as a particle rotating around X2 and 
have a nonzero charge radius. However, for the mass degenerate case, X^ and X2 should be treated 
with equal foot and rotate around the center with the same radius. As a result, for each orbit the net 
charge is zero and the charge radius is zero. 

Similarly, we can integrate out Xi and X2 to generate the effective operator coupling S to two 
photons. To match the coefficient in Eq. ([7D, we have 

M 12 [mj,^ m\j ' 



In the limit rnj^^/fi^ <^ m^^/^2) we have 



Using the values of M/Xs in Fig. HI we anticipate at least the charged particle X2 to have a mass 
below 1 TeV. This charged particle X2 can decay into one lepton plus one neutrino, for example via 
the higher dimensional operator X2HLeR. 



^Operator A^Xixe^ does not preserve this symmetry, but this operator could have a very small coefficient and is 
irrelevant to this calculation. 
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3.3 Dark Matter Relic Abundance 

In our DeDM model, we have two stable particles in our spectrum: x (pi. In our previous 
analysis, we have assumed that the majority of dark matter in our universe is mainly composed of x- 
To justify our assumption, it is important to study the thermal history of x ^-iid (pi. In this section, 
we demonstrate that our setup contains enough ingredients to induce a right relic abundance for Xi 
thus it could be the dominant part of the dark matter in our current Universe. 

The thermal relic abundance of (pi is controlled by the parameter Ai in Eq. (I15p . which is similar 
to the "Higgs portal" dark matter models [39|[50]. For m^^ < m/^, the main annihilation cross section 

is EH 



f^Mn) - 77712 ,2N2 . --.-2t.2/„.. X ^777 ' i-^") 



where vew = 246 GeV is the electroweak vacuum expectation value. The function Tfi{m) is the width 
of a Higgs boson in the SM with a mass at m. For Ai = 1, m^^ = 100 GeV and nih = 125 GeV, we 
have avr{(pi) ~ 581 pb and ^^^h? ~ 1.4 x 10^^ x Q.-Qyih? . Thus the relic abundance of <pi can be 
naturally small. 

To satisfy the dark matter relic abundance, a non-trivial thermal history of x is needed. This 
is because a large annihilation cross section in Eq. (jlOp is needed to explain the INTEGRAL data. 
The thermal relic abundance of x is very small compared to the required dark matter energy density. 
Noticing that the last operator in Eq. (jl5p can introduce the decay channel, — )■ xe"*", the late 
decay of thermally abundant X^ particles can generate enough x, and therefore explain why x could 
be the majority of dark matter. 

We first calculate the thermal relic abundance of the charged particle X^ before it decays into 
X and a positron/electron. There are two classes of annihilation channels for X^ . The first class 
has a photon or Z boson exchanging in the s-channel with final states as a pair of the SM fermions, 
W~ gauge bosons, and h Z. The second class includes the t-channel diagrams, interfering with 
seagull diagrams. Assuming that the mass X^ is far above the SM particle masses and neglecting the 
SM particle masses, we have 

.„„WXr^7T) = ...WXr^ZZ) = ^^i^, (21) 
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aVr{XtX^ ^W+W~) = -vl, (23) 

1536 vr Cy^ m^_^ 



4 / 

e nc 



az;,(X+Xr ^ //) = 2 ' 2 " " '^^r, (24) 



96 vr cl 



Here, g-^ is the electric charge of the SM fermion; g-^ iOv) is the axi-vector (vector) couphngs of the 
Z to the SM fermion up to the electric coupling e; ni = 3 for quarks and 1 for leptons. To derive 
the above formulas, we have only included the leading terms in Vr for each equation. If the charged 
particle could occupy the total energy of dark matter in our universe, J^dm^^ = 0.11, the required mass 
is calculated to be ~ 750 GeV. To derive this mass, we have found that the p-wave suppressed 
annihilation cross section or the terms at 0{v'^) is subdominant compared to the total annihilation 
cross section. 

For a lifetime of not too long in the cosmological time scale, we should anticipate that X^ 
has already decayed into its daughter particle and the final dark matter in our current universe is 
composed of x- On the other hand, the lifetime of X^ can not be too short. Otherwise, the produced 
X particles from Xf^ decays in the early universe can easily annihilate away and do not provide enough 
dark matter energy density. To calculate the thermal history of the x field, one needs to solve for the 
following coupled Boltzmann equations between x and X^ 

+ 3Hnx^ = -{(Tv)x^ (n^i - n^^) - "-Xi , (25) 

diTi/ 

^ + SHn^ = - {av)^ - n^^^) + nx. Fx, . (26) 

Here, in the radiation dominated era, H = (87r/9/3Mpi)i/^ t = 1/{2H), p{T) = 5* vr^ r'^/30, and 
n^^{T) = gfj(mjT/27r)'^/^e~"^'/"^, where = 86.25 is the number of degrees of relativistic freedom and 
gi^- = 4 and gxi =2. It is convenient to rescale the number density by the entropy and to define the 
quantity Yi = rii/s with s = 27r^(7*T^/45. The coupled equations become 

where x = fn^/T and dx/dt = Hx. The final x relic abundance is given by = Px/ Pc^ where 
Pc = 3HqMp^/8-it = 1.0539 x lO^^/i^ GeVcm^''^ is the critical density corresponding to a flat universe 
and Px = mxSoYx{oo) with sq = 2889.2 cm~^ being the entropy today. 

At the temperature region with 0(20) < x < 1000, the decaying terms in Eq. (j27|) and (|28|) are not 
important. The number densities of Xi and x reach their separate freeze-out values. Since the cross 
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Figure 5: The comoving number density as a function of the temperature. Here, we choose 
rriy. = 135 GeV and mxi = 3.8 TeV. The annihilation cross section of Xi, approximately 0.03 pb, is 
determined by its interactions with electroweak gauge bosons. The solid lines are for {avr)x = 10 pb 
and Txi = 10 s; the dotdashed lines are for {crvr)x = 10 pb and txi = 100 s; the dotted and red line 
is for {(jVr)x = 100 pb and tx-^ = 10 s. The relic abundance of x satisfies the observed dark matter 
energy density, ^^.^'^ ~ 0.11- 



section of x is much larger than Xi, the freeze-out number density for x is much below the one of Xi. 
At a later time, only the last terms in Eq. (I27p and (I28|) become important. One can easily show that 
the quantity Yx^ + 1^ is a conserved number. As a result, the final number density of x should just 
match to the number density of Xi at 0(20). So, approximately we have the relic abundance of x &s 

l^./i^^O.ll X X (29) 

^ 135 GeV 3.8 TeV ^ ' 

So, the charged particle Xi is predicted to be 3.8 TeV and the other charged particle X2 should be 

below around 1 TeV to explain the ratio of INTEGRAL and Fermi-LAT cross sections in Eq. (|18p . 

We solve the coupled equations in Eq. (|28p numerically and show both comoving number densities 

of X s-iid ^1 ill Fig- El We find that if the mass of Xi is 3.8 TeV, it generates the relic abundance 

for X which satisfies the total dark matter energy density, ^^^^^ ~ 0.11. In the blue solid and the red 

solid lines, for {aVr)x = 10 pb[f| and txi = 10 s we show the evolutions of the Xi and x comoving 

number densities as a function of temperature. As can be seen from Fig. [5] and at x ~ 20, both Xi 

and X have reached ordinary relic abundances according to their respective annihilation cross sections. 



^The annihilation cross section of x is not necessarily related to its annihilation cross section at the current time. 
This is because its main production here is from the heavy particle Xi decay, and it has a relativistic velocity and hence 
a smaller cross section. 
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At X ~ 10^ — 10^, Xi starts to decay and its number density drops rapidly. Meanwhile, the stable x 
particle number density increases and reaches a plateau at around 10''. The final number density of 
the X field is found to be independent on the lifetime tx-^ , as long as the decay happens late enough 
so that the annihilation of x is not important any more. The actual time for x to reach its eventual 
number density is proportional to jTx-y- 




To satisfy the dark matter relic abundance, the lifetime of the charged particle X\ can be < 
100 s . For such a late decayed particle, we need to worry about its modification on the Big Bang 
nucleosynthesis (BBN) history. Since the main decaying product of X\ is into leptons plus the stable 
X, the BBN constrains are fairly weak. From Ref. [52], the ^Li/^Li ratio constrains the lifetime of X\ to 
be < 10^ s for ^Xx^ ^ 0.1 if it would have not decayed. As pointed in Ref. |53U54j . the long-lived 
charged particle, with a lifetime tx^ > 10^ s, can form a bound state with nuclei and enhance the ^Li 
production. The parameter space in our model can indeed satisfy the BBN constraints. 

4 Discussion and conclusions 

The charged particle Xi in our model behaves as a heavy stable charged particle (HSCP) at colliders. 
The current searches from CMS at ^/s = 7 TeV and 5.0 fb^^ have set a lower limit on the mass of Xi 
to be 223 GeV at 95% C.L. [55]. For the predicted mass of Xi around 3.8 TeV, the existing studies 
have shown that the 14 TeV LHC with 100 fb~^ can reach the HSCP up to a mass around 1 TeV. 
So, unlikely the stable charged particle can be discovered at the 14 TeV LHC. However, for the other 
charged particle X2 its mass should be below 1 TeV and could be a long-lived particle or decay into 
SM particles, for instance — )• e+fg. The parameter space of the X2 particle will be well explored 
at the LHC 14 TeV running. 

One feature of our model is directly using photon as a mediator to link the dark matter sector 
to positron/electron. Unfortunately, other than searching for the charged particles responsible for 
the charge radius operator, in the near future there is no additional observable dark matter direct or 
indirect signatures for the x field, which has interactions with SM particles suppressed by the TeV 
scale cutoff of the effective operators. The minor component of dark matter, 0i, may have detectable 
effects. However, that highly relies on the parameters one chooses, thus we do not pursue that in 
detail here. Another ingredient that we utilize is the s-channel resonance particle to increase the 
annihilation cross section. We want to stress that this option is not a unique one and is introduced 
just for convenience. One can also introduce a light mediator in the i-channel plus the Sommerfeld 
enhancement to achieve the same goal [48| . 

In summary, we have constructed a realistic model to have the same dark matter particle responsi- 
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ble for both the INTEGRAL 511 keV and Fermi-LAT 135 GeV Unes. Through an s-channel resonance, 
the dark matter particles annihilate into a complex scalar, which couples to photon via a charge-radius 
operator. For a few MeV mass splitting between the real and imaginary parts of the complex scalar, 
two pairs of electron and positron are the main visible particles from dark matter annihilation. We 
have worked out the parameter space and have found that both the large cross section required for 
INTEGRAL and the small cross section for Fermi-LAT can be simultaneously accommodated in our 
model. The thermal relic abundance of dark matter is achieved by the late decay of a charged par- 
ticle, which also generates the charge-radius operator. The other charged particle responsible for the 
charge-radius operator is predicted to have a mass below 1 TeV. The 14 TeV LHC will concretely test 
the scenario presented in this paper. 
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